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Abstract 

The hodographs of magnetization of nonlinear NMR are investigated in the 
conditions of resonance on the unshifted frequency. It is shown that, depend- 
ing on the value of amplitude of the variable field and value of frequency shift, 
topologically different hodographs separated from each other by separatrix 
are obtained. It is shown that the set of hodograph points, being obtained 
by the stroboscopic method, is chaotic and the change of ^-component of 
magnetization has the form of solitons chaotically changing the sign. 
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1. Introduction 

Investigations of nonlinear magnetic resonance in magnetically ordered 
materials have been carried out for a long time [1], however the interest 
to them is still increasing due to discovery of new physical phenomena in 
magnetics [2-9]. On the other hand the interest to chaotic Hamiltonian 
systems does not weaken [10-12]. That is why the study of any questions 
on the butt of these two fields (our work concerns to) in our opinion is 
quite topical. In particular most of spin systems considered as promising 
candidates for applications in quantum computation are nonlinear [13-17]. 
Nonlinearity is a principle problem for controlled states. Note that for a linear 
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system in a constant magnetic field, magnetization vector can be deflected 
by applying a variable radio frequency field. However for nonlinear systems, 
situation is completly different. 

In magnetically ordered crystals at quite, low temperatures a shift of 
frequency precession of nuclear magnetization m proportional to the longi- 
tudinal component m z (dynamical shift) origins. In the NMR conditions 
appearance of dynamic shift leads to essentially nonlinear character of the 
motion m. For example, in [18] it was shown, that fulfilling the conditions 
uj\ = Up, where oj p is equilibrium value of dynamic shift, u)\ = rj^Hi, Hi 
is an amplitude of linearly polarized variable field, applied in the transverse 
plane, r\ is gain factor, 7 is a gyromagnetic ratio, the motion of m vector in 
the rotational coordinate system becomes aperiodic (period T — > 00). Tra- 
jectory, along which the end of m vector moves (hodograph) in this case is a 
separatrix topologically separating differing trajectories. It will be shown be- 
low that points put according to hodograph being close to separatrix form a 
stochastic set. Simple, but effective method of finding the stochastic motion 
(stroboscopic method) is offered. 

Peculiarity of dynamical systems is extreme sensitivity of motion near 
separatrix, even with respect to the slightly change of initial conditions or 
adiabatic perturbation. Purpose of our paper is demonstration at the fact 
that control of the state of nonlinear system is not possible near separatrix. 
In order to proof this statement, we shall consider realistic physical system. 
Namely, we shall study motion of nuclear magnetization in the presence of 
large dynamical shift of frequency. 

2. Equations of motion of magnetization at nonlinear NMR 

Let us direct the axis z along the total field H = H + H; oc , where H and 
H; oc are respectively external and internal magnetostatic fields on the nuclei. 
Let us assume that the variable fields also effect on the nuclear spin-system 

Hi (t) = Hi cos cut, Hi±H (1) 

and 

hi(i) = hicosi/t, h || H (2) 

Directing Hi along the axis x, the system of nonlinear differential equations 
describing the motion of nuclear magnetization with the account of dynamical 



2 



shift of frequency precession in the rotating system of coordinates around the 
axis z with the frequency can be written as follow^ 



x — (A — u) p z)y + v\y cos vt, 



y = —(A — uo p z)x + 2 UlZ ~ VlX cos 
z = —-ooiy, v\ = 'jh, A = u — u , ujq = jH. 



(3) 



Here we introduced dimensionless components of magnetization a = m a /m, 
where m = |m|, a = x,y,z. Our observation neglects relaxation effects 
and magnetostatic inhomogeneity of the local field is valid in the case of 
fields with duration r <^ Ti, T 2 , 2n/Q, where 7\ and T 2 are the times of 
longitudinal and transverse relaxation, Q is inhomogeneous width of the line. 
The system (3) is investigated in the conditions of weakness of longitudinal 
fields (z/x uj p ). 

We investigate an unperturbed system being obtained from (3) at v\ = 0. 
It is easy to see that unperturbed system has two integrals of motion: 



The first of these integrals (4) corresponds to conservation of the value 
of full magnetic moment, the second (5) corresponds to the conservation of 
full magnetic energy. Here we assume that the conditions of resonance are 
fulfilled on the unshifted frequency A = 0, and x(0) = y(0) = 0, z(0) = 1 
are taken as initial conditions. These relations give definite surfaces: (4) 
single sphere and (5) parabolic cylinder with generatrix along the axis y. 
Hodographs of vector m are closed curves obtained by crossing of these sur- 
faces. 

When uji ^ uj p , hodographs are divided into two types: 1) U\ < u p - the 
section consists of two contours, which are symmetrically located on different 
sides of the plane (x,y); 2) lj 1 > oj p - the section consists of one contour 



1 Investigation of Bose-Hubbard model is also reduced to consideration of such type of 
equations [19] 



x 2 + y 2 + z 2 = 1 

UJpZ 2 — LO\X = Up. 



(4) 
(5) 
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twice crossing the plane (x,y). When cu 1 = oj p the trajectory looks like 
symmetrical spatial eight with self-crossing on the at x — — 1 (Fig.l). 

One-parametric (parameter uji) family of eights creates the separatrix of 
unperturbed motion, which divides the space into domains according to the 
type of their trajectories. 




Figure 1: Family of hodographs of vector r = m/|mo| in the conditions of resonance on 
the unshiftcd frequency: 1) u>i < u> p . Two closed curves, located symmetrically relatively 
to the plane (x,y). 2) u)\ > u) p . The curve looks like the spatial eight, which crosses the 
plane (x, y) twice. 3) uo\ = uj p . The spatial eight with self-crossing on the axisx at x = — 1 
acts as separatrix. 



Using the integrals of motion (4) and (5) unperturbed (v\ = 0) system of 
equation (3) can be reduced to 



z = -2 U p 



-{l- Z 2 )[l-^- Z 



Solution of equation (6) can be written as: 

dn\t'\uji/uj^[ uj\ < ijj. 



1 cn[t ;ui/u p \ uji > uj p 



1/2 



(6) 



(7) 



t' — —-uj p t, 

where cn and dn are elliptical functions of Jacob, cosine and delta of ampli- 
tude. 
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For the periods of motion at closed trajectories in the case of U\ < oj p , 
we can obtain relatively 

T_= f-Vf-Y ^<"v (8) 

and 



Up/ \iO p 



where K(k) is full elliptical integral of the first order. Hence for the period 
of motion close to the separatrix lo\ ~ u p we obtain 

T c = T+ « 2T « — Zw . 4u;i , . (10) 
^i |^i ~~ 

With the approximation to the separatrix the period of motion logarith- 
mically diverges. The right part of equation (3) at the point (-1,0,0), near 
the knot of eight, becomes zero. It means, the motion moderates near this 
point. It is evident that such motion takes place close to separatrix hodo- 
graph (Fig.l). Thus, close to the separatrix the motion is uneven. Near the 
point (-1,0,0) it moderates, and the remaining part of hodograph is quite 
rapidly overcome, so that full time of one rotation is determined by the ex- 
pression (10). Time dependence of z component of magnetization will have 
the form of instanton [20]. The uj\ = u p determines the bifurcation value of 
amplitude of the variable field. 

Note, that at A = condition of resonance is fulfilled not at the start of 
motion (point (0,0,1)), but at reaching special point (-1,0,0). It is clear that 
the peculiarities of the motion near the special point are caused by fulfilling 
the resonance condition in this point. 

Now we discuss the case of resonance on the shifted frequency A — uj v = 0. 
In this case the conditions of resonance are fulfilled at the start of motion in 
the point (0,0,1). In NMR experiments in the materials with big dynamic 
shift, resonance condition is selected, namely, in such a way. In this case the 
integrals of motion have the form: 

x 2 + y 2 + z 2 = l (11) 

UJ p (l - Z f = LU X X (12) 

It is clear that the crossing contours of two surfaces (11) and (12) are 
topologically equitype closed curves (Fig. 2), that is why from dynamical 
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stochasticity point of view it is not interesting. Below we consider only the 
case of resonance on the unshifted frequency (A = 0). 




Figure 2: The family of hodographs of vector r = m/|m | in the conditions of resonance 
on the shifted frequency A — oj. p = 0. With the growth of w\{(Jp < lo^ < . . . the 
hodographs increase in dimensions and in the limit lo\ « lo p reach the point (-1,00). 



3. Analysis of Results of Computer Modeling 

For analysis of the above mentioned discussions and investigation of the 
influence of perturbation on motion of nuclear magnetization the initial equa- 
tions (1) describing spin-system dynamics were modeled. Modeling was made 
in the medium of Simulink of mathematical package MATLAB. During the 
modeling we used stroboscopic method, according to which observation of 
vector motion of magnetization was conducted not continuously, but during 
separate, short periodically following one another time intervals. 

Fig. 3 shows hodographs distant from separatrix {uj 1 ^ uj p ). The points 
being set using stroboscopic method are periodically located along the hodo- 
graph. Approaching the separatrix (ui — > uo p ) they become chaotic (Fig.4). 
In the case of uj\ ~ uo p and v\ — chaotic points cover only the upper part of 
separatrix eight, but after switching on perturbation [y-y = \Q~ A Hz) the mo- 
tion acts along the whole separatrix one. Stroboscopic points of hodograph 
are chaotic in this case too. Consequently, a role of slowly changing pertur- 
bation is reduced to regular variation of the place of the point of self-crossing 
separatrix eight. 
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Figure 3: Stroboscopic points of hodographs, which are far from separatrix (wi ^ 
Up), a) wi = 800Hz, uj p = lOOOiJz, v x = 0; b) w\ = lOOOTJz, uj p = 800iJz, v x = 0. In 
both cases the set of stroboscopic points is located periodically. 

Thus, stochastic set of points origins close to the singular point (-1,0,0), 
overcoming of which the motion becomes unpredictable: all four directions, 
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Figure 4: Stroboscopic points of separatrix hodograph: a) u\ = u p = 1000-ffz, v\ = 0. 
Chaotic points located on the upper part of hodograph. b) u>\ — ui p = 1000Hz, v\ = 
10~ 4 Hz. Chaotic points located along the whole separatrix hodograph. 



crossing the special point become equally probable. Notice that set of points 
turned out to be stochastic though external influence on the system is peri- 



odical. 

As a result of numerical integration of (3) for z(t) at the motion close to 
the separatrix the solution in the form of periodical succession of instantons 
is obtained (Fig. 5). In the case of lo\ — u p — > —0 all the instantons are 
positive, while at lo\ — uj v — > +0 in the periodical succession alternate change 
of instanton sign takes place. 

The same result can be obtained also by means of direct set of figures of 
elliptical functions of Jacob (7). 




Figure 5: Change of z component of nuclear magnetization during simulation for separatrix 
close to trajectory a) ui p > wi (u)\ = 1000Hz, ui p = 1001Hz, v\ =0); b) ui p < ui\ (wi = 
1001Hz, Up = lOOOiJz, v x = 0). 

Fig. 6 shows the fragment of the picture of changing the sign of m z during 
simulation already at the switched periodical perturbation {y\ = 10~ A Hz). 
Disordered change of the sign of nuclear magnetization testifies an origin of 
stochasticity at the junction over branching point. 

Thus, hodograph of vector of nuclear magnetization at nonlinear NMR 
on the unshifted frequency in the conditions lo\ ~ uo p possess the properties 
of separatrix. This is the basis for appearance of the set of chaotic points in 
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Figure 6: Change of z component of vector of nuclear magnetization during simulation for 
separatrix trajectory (ui p — u>\ = lOOOiJz, v\ = 10~ A Hz). 

the stroboscopic picture of separatrix hodograph (Fig. 3), and also for chaotic 
change of instantons signs of z component of magnetization (Fig. 6). 
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